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We theoretically show that intriguing features of coherent many-body physics can be observed in 
electron transport through a quantum dot (QD). We first derive a master equation based framework 
for electron transport in the Coulomb-blockade regime which includes hyperfine (HF) interaction 
with the nuclear spin ensemble in the QD. This general tool is then used to study the leakage current 
through a single QD in a transport setting. We find that, for an initially polarized nuclear system, 
the proposed setup leads to a strong current peak, in close analogy with superradiant emission of 
photons from atomic ensembles. This effect could be observed with realistic experimental parameters 
and would provide clear evidence of coherent HF dynamics of nuclear spin ensembles in QDs. 



I. INTRODUCTION 

Quantum coherence is at the very heart of man y in- 
triguing phenomena in today's nanostructures^l. For 
example, it is the essential ingredient to the understand- 
ing of the famous Aharonov-Bohm like interference os- 
cillations of the conductance of metallic rings^ or the 
well- kno wn conductance steps in quasi one-dimensional 



wirea^l. In particular, nonequilibrium electronic trans- 
port has emerged as a versatile tool to gain deep insights 
into the cohe ren t quantum properties of mesoscopic solid- 
state devices^'. Here, with the prospect of spintron- 
ics and applications in quantum computing, a great deal 
of research has been directed towards the interplay and 
feedback mechanisms between electron and nuclear spins 
in gate-based semiconductor quantum dots^H^. Current 
fluctuations have been assigned to the random dynam- 
ics of the ambient nuclear spins^ and/or hysteresis ef- 
fects due to dynamic nuclear polarizatio n 15 " 18 !. Spin- 
flip mediated transport, realized in few-electron quantum 
dots in the so-called spin-blockade regime^, has been 
shown to exhibit long time scale oscillations and bista- 
bility as a result of a buildup and relaxation of nuclear 
polarizatior i 15 * 16 *. The nuclear spins are known to act col- 
lectively on the electron spin via hyperfine interaction. In 
principle, this opens up an exciting testbed for the obser- 
vation of collective effects which play a remarkable role 
in a wide range of many-body physics^^l. 

In Quantum Optics, the concept of superradi- 
ance, describing the cooperative emission of photons, 
is a parad igm example for a cooperative quantum 
effeclfjJMH] Here, initially excited atoms emit photons 
collectively as a result of the buildup and reinforcement of 
strong interatomic correlations. Its most prominent fea- 
ture is an emission intensity burst in which the system ra- 
diates much faster than an otherwise identical system of 
independent emitters. This phenomenon is of fundamen- 
tal importance in quantum optics and has been studied 
extensively since its first prediction by Dickc in 1954 23 . 
Yet, in its original form the observation of optical super- 
radiance has turned out to be difficult due to dephasing 
dipole-dipole van der Waals interactions, which suppress 
a coherence buildup in atomic ensembles. 



This paper is built upon analogies between mesoscopic 
solid-state physics and Quantum Optics: the nuclear 
spins surrounding a QD arc identified with an atomic 
ensemble, individual nuclear spins corresponding to the 
internal levels of a single atom and the electrons are as- 
sociated with photons. Despite some fundamental dif- 
ferences - for example, electrons are fermions, whereas 
photons are bosonic particles - this analogy stimulates 
conjectures about the potential occurence of related phe- 
nomena in these two fields of physics. Led by this line 
of thought, we will address the question if superradiant 
behaviour might also be observed in a solid-state envi- 
ronment where the role of photons is played by electrons. 
To this end, we analyze a gate-based semiconductor QD 
in the Coulomb blockade regime, obtaining two main 
results, of both experimental and theoretical relevance: 
First, in analogy to superradiant emission of photons, we 
show how to observe superradiant emission of electrons in 
a transport setting through a QD. We demonstrate that 
the proposed setup, when tuned into the spin-blockade 
regime, carries clear fingerprints of cooperative emission, 
with no van der Waals dephasing mechanism on relevant 
timescales. The spin-blockade is lifted by the HF cou- 
pling which becomes increasingly more efficient as corre- 
lations among the nuclear spins build up. This markedly 
enhances the spin-flip rate and hence the leakage current 
running through the QD. Second, we develop a general 
theoretical master equation framework that describes the 
nuclear spin mediated transport through a single QD. 
Apart from the collective effects due to the HF interac- 
tion, the electronic tunneling current is shown to depend 
on the internal state of the ambient nuclear spins through 
the effective magnetic field (Overhauser field) produced 
by the hyperfine interaction. 

The paper is structured as follows: In Sec.|ITJ we high- 
light our key findings and provide an intuitive picture of 
our basic ideas, allowing the reader to grasp our main 
results on a qualitative level. By defining the underlying 
Hamiltonian, Sec. |III| then describes the system in a more 
rigorous fashion. Next, we present the first main result 



of this paper in Sec. IV a general master equation for 
electron transport through a single QD which is coher- 
ently enhanced by the HF interaction with the ambient 
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FIG. 1. (color online). Schematic illustration of the trans- 
port system: An electrically defined QD is tunnel-coupled to 
two electron reservoirs, the left and right lead respectively. A 
bias voltage eV = flL — PR is applied between the two leads 
in order to induce a current through the QD. An external 
magnetic field is used to tune the system into the sequential- 
tunneling regime and the QD effectively acts as an spin- filter. 
The resulting spin-blockade can be lifted by the HF interac- 
tion between the QD electron and the nuclear spins in the 
surrounding host environment. 



nuclear spins in the QD. It features both collective effects 
and feedback mechanisms between the electronic and the 
nuclear subsystem of the QD. Based on this theoretical 
framework, Sec. [V] puts forward the second main result, 
namely the observation of superradiant behavior in the 
leakage current through a QD. The qualitative explana- 
tions provided in Sec. [Il] should allow to read this part 
independently of the derivation given in Sec. |IV| Sec. |VI| 
backs up our analytical predictions with numerical sim- 
ulations. When starting from an initially polarized nu- 
clear spin ensemble, the leakage current through the QD 
is shown to exhibit a strong peak whose relative height 
scales linearly with the number of nuclear spins, which 
we identify as the characteristic feature of superradiant 
behaviour. In Sec. |VII] we draw conclusions and give an 
outlook on future directions of research. 



II. MAIN RESULTS 

In this section we provide an intuitive exposition of our 
key ideas and summarize our main findings. We study 
a single electrically-defined QD in the Coulomb-blockade 
regime which is attached to two leads, as schematically 
depicted in Fig. [T] Formally, the Hamiltonian for the 
total system is given by 



H — Hz + Hb + Ht + H) 



HF- 



(1) 



Here, Hz describes the Zeeman splitting of the electronic 
levels inside the QD in presence of an external magnetic 



field. Hb refers to two independent reservoirs of non- 
interacting electrons, the left and right lead respectively. 
The coupling between these and the QD is described in 
terms of a tunneling Hamiltonian Ht and i?HF models 
the collective interaction between an electron confined 
inside the QD and an ensemble of N proximal nuclear 
spins surrounding the QD. Note that the specific form of 
H will be given later on in Sec. |III| 

Our analysis is built upon a Quantum master equa- 
tion approach, a technique originally rooted in the field 
of Quantum Optics. By tracing out the unobserved de- 
grees of freedom of the leads we derive an effective equa- 
tion of motion for the density matrix of the QD system 
ps - describing the electron spin inside the QD as well as 
the nuclear spin ensemble - irreversibly coupled to source 
and drain electron reservoirs. In addition to the standard 
assumptions of a weak system-reservoir coupling (Born 
approximation), a flat reservoir spectral density, and a 
short reservoir correlation time (Markov approximation) , 
we demand the hyperfine flip-flops to be strongly detuned 
with respect to the effective magnetic field seen by the 
electron throughout the dynamics. Under these condi- 
tions, the central master equation can be written as 



p s (t) = -i[H z +H nF ,p s (t)] 



(2) 



daPS (t)4r~ -{dld a ,p S (*)} 



+ M*) 4ps(*K-5{<k4>As(*)} 

<r=t,4. L 

where the tunneling rates a a (t) and /3 CT (t) describe dis- 
sipative processes by which an electron of spin a tun- 
nels from one of the leads into or out of the QD, respec- 
tively. Here, the fermionic operator d\ creates an electron 
of spin a inside the QD. While a detailed derivation of 
Eqn.Q along with the precise form of the tunneling rates 
is presented in Sec. IV here we focus on a qualitative dis- 



cussion of the theoretical and experimental implications 
thereof. 

Our central master equation exhibits two core features: 
First, dissipation only acts on the electronic subsystem 
with rates a a (t) and j3 a (t) that depend dynamically on 
the state of the nuclear subsystem. This non-linear be- 
havior potentially results in hysteretic behavior and feed- 
back mechanisms bet ween the tw o subsystems as already 
suggested theoreticall y^ 11 * 14 * 20 * 21 ! and observed in experi- 
ments in the context of doubl e QDs in the Pauli-blockade 
regime; see e.g. RefsJ^nSESl Second, the collective na- 
ture of the HF interaction Hup allows for the observation 
of coherent many-body effects. 

The effect of the hyperfine interaction between an elec- 
tron inside the QD and the ambient nuclear spin ensem- 
ble is two-fold giving rise to the two main results out- 
lined above: First, the nuclear spins provide an effec- 
tive magnetic field for the electron spin, the Overhauser 
field, whose strength is proportional to the polarization 
of the nuclear spin ensemble. Thus, a changing nuclear 
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polarization can either dynamically tune or detune the 
position of the electron levels inside the QD. This, in 
turn, can have a marked effect on the transport proper- 
ties of the QD as they crucially depend on the position of 
these resonances with respect to the chemical potentials 
of the leads. In our model, this effect is directly cap- 
tured by the tunneling rates dynamically depending on 
the state of the nuclei. Second, to show that this system 
supports the observation of intriguing, purely collective 
effects we refer to the following example: Consider a set- 
ting in which the bias voltage and an external magnetic 
field are tuned such that only one of the two electronic 
spin-components, say the level lies inside the trans- 
port window. In this spin-blockade regime the electrons 
tunneling into the right lead are spin-polarized, i.e., the 
QD acts as an spin filter^EH. If the HF coupling is suffi- 
ciently small compared to the external Zeeman splitting, 
the electron is predominantly in its |4-) spin state allowing 
to adiabatically eliminate the electronic QD coordinates. 
In this way we obtain an effective equation of motion for 
the nuclear density operator p only. It reads 



A~pA* 



1 



{A+A-,p} 



+id [A+A ,n] +i| [A z ,p] , 



(3) 



where A^ — J^iLi 9iC^ with fi = +,—,z are collective 
nuclear spin operators, composed of all N individual nu- 
clear spin operators erf , with gi being proportional to 
the probability of the electron being at the location of 
the nucleus of site i. Again, we will highlight the core 
implications of Eqn.(|3| and for a full derivation thereof, 
including the definition of the effective rates c r and 
we refer to Sec. [V] Most notably, Eqn.(|3| closely resem- 
bles the superradiance master equation which has been 
discussed extensively in the context of atomic physics^ 
and therefore similar effects might be expected. 

Superradiance is known as a macroscopic collective 
phenomenon which generalizes spontaneous emission 
from a single emitter to a many-body system of N 
atoms^. Starting from a fully polarized initial state the 
system evolves within a totally symmetric subspace un- 
der permutation and experiences a strong correlation 
build-up. As a consequence, the emission intensity is 
not of the usual exponentially decaying form, but con- 
versely features a sudden peak occuring on a very rapid 
timescale ~ 1 /N with a maximum ~ N 2 . 

In this paper, we show that the same type of coopera- 
tive emission can occur from an ensemble of nuclear spins 
surrounding an electrically-defined QD, a phenomenon 
we term as electronic superradiance: The spin-blockade 
can be lifted by the HF interaction as the nuclei pump ex- 
citations into the electron. When starting from a highly 
polarized, weakly correlated nuclear state, due to the col- 
lective nature of the HF interaction, this process becomes 
increasingly more efficient as correlations among the nu- 
clei build up, directly giving rise to an increased leakage 
current. Therefore, the current is collectively enhanced 
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FIG. 2. (color online). Normalized leakage current through 
a QD in the spin-blockade regime for N nuclear spins, initial 
nuclear polarization p and external Zeeman splitting cjo in 
units of the total HF coupling constant Ahf ~ 100 fieV, sum- 
marized as (N,p, ljo/Ahf)- For homogeneous HF coupling 
the dynamics can be solved exactly (black dotted line) . Com- 
pared to this idealized benchmark, the effects are reduced for 
realistic inhomogeneous HF coupling, but still present: The 
relative peak height becomes more pronounced for smaller 
detuning uiq or higher polarization p (solid red line compared 
to the blue dashed and green dash-dotted line, respectively). 
Even under realistic conditions, the relative peak height is 
found to scale linearly with N, corresponding to a strong en- 
hancement for typically N w 10 5 — 10 6 . 



by the electron's HF interaction with the ambient nuclear 
spin ensemble giving rise to the pure many-body effect 
of electronic superradiance. 

Compared to its conventional atomic counterpart, our 
system incorporates two major differences: First, our 
setup describes superradiant behaviour from a single 
emitter, since in the strong Coulomb-blockade regime the 
electrons are emitted antibunched. As described above, 
the superradiant character is due to the nuclear spins 
acting collectively on the electron spin leading to an in- 
creased leakage current on timescales longer than single 
electron tunneling events. The second crucial difference 
is the inhomogeneous nature (<7, ^ const.) of the collec- 
tive operators A* 1 . Accordingly, the collective spin is not 
conserved, leading to dephasing between the nuclei which 
in principle could prevent the observation of superradi- 
ant behavior. However, as exemplified in Fig. [2j we show 
that under realistic conditions - taking into account a 
finite initial polarization of nuclear spins p and dephas- 
ing processes due to the inhomogeneous nature of the HF 
coupling - the leakage current through the QD still ex- 
hibits the characteristic peak whose relative height will 
be shown to scale linealry with the number of nuclear 
spins. Even though the effect is reduced compared to the 
ideal atomic case, for an experimentally realistic number 
of nuclei N s» 10 5 — 10 6 a strong increase is still predicted. 
The experimental key signature of this effect, the relative 
peak height of the leakage current, can be varied by ei- 



4 



ther tuning the external Zeeman splitting or the initial 
polarization of the nuclear spins. 

In the remainder of the paper, Eqn.Q and Eqn.([3| 
will be derived from first principles; in particular, the un- 
derlying assumptions and approximations will be listed. 
Based on this general theoretical framework, more re- 
sults along with detailed discussions will be presented. 
For both the idealized case of homogeneous HF coupling 
- in which an exact solution is feasible even for relatively 
large N - and the more realistic inhomogeneous case, 
further numerical simulations will prove the existence of 
a strong superradiant peaking in the leakage current of 
single QD in the spin-blockade regime. 



III. THE SYSTEM 

This section gives an in-depth description of the Hamil- 
tonian under study, formally introduced in Eqn.Q. The 
system we consider consists of a single electrically-defined 
QD in a transport setting as schematically depicted in 
Fig. [I] Due to strong confinement only a single or- 
bital level is relevant. Moreover, the QD is assumed 
to be in the strong Coulomb-blockade regime so that 
at maximum one electron resides inside the QD. There- 
fore, the effective Hilbert-space of the QD electron is 
span , ||) , |0)} where the lowest energy states for an 
additional electron in the QD with spin a =t, i arc split 
by an external magnetic field. The Hamiltonian for the 
total system is given in Eqn. ([!]). 

Here, the first term H z = J2 a e <jd\d<j describes the 
electronic levels of the QD. The Zeeman splitting be- 
tween the two spin components is ujq — — (we set 
H = 1) and the QD electron operators are d\ = \a) (0|, 
describing transitions from the state |0) with no electron 
inside the QD to a state |<r) with one electron of spin a 
inside the QD. 

Electron transport through the QD is induced by at- 
taching the QD to two electron leads (labeled as L and R) 
which are in thermal equilibrium at chemical potentials 
and fin, respectively. The leads themselves constitute 
reservoirs of non-interacting electrons 



tonian 



e otk c a kcr C cika, 

Oi.k.a 



(4) 



where c? aka (c a k<j) creates (annihilates) an electron in 
lead a = L, R with wavevector k and spin a. The 
operators c^ fc(T (c a ka) fulfill the usual Fermi commuta- 
tion relations: {c\ ktT , 4,/ fcv ,} = {c Qfccr , cv fc ^} = 

and { c akcr' c a'k>a>} = &a,ot' &k,k' ■ The effect of the 
Coulomb interaction in the leads can be taken into ac- 
count by renormalized effective quasi-particle masses. A 
positive source-drain voltage eV = /il — fiR leads to a 
dominant tunneling of electrons from left to right. Micro- 
scopically, the coupling of the QD system to the electron 
reservoirs is described in terms of the tunneling Hamil- 



(5) 



with the tunnel matrix element T£ a J. specifying the trans- 
fer coupling between the lead a = L 7 R and the system. 
There is no direct coupling between the leads and elec- 
tron transfer is only possible by charging and discharging 
the QD. 

The cooperative effects are based on the the collective 
hyperfine interaction of the electronic spin of the QD 
with N initially polarized nuclear spins in the host envi- 
ronment of the QEBH. It is dominated by the isotropic 
contact terrrPfil given by 



H 



HF 



+ A-S+) +gA z S z 



(6) 



Here S* 1 and = J2iLi9i a i wrtri M = +j — i z de- 
note electron and collective nuclear spin operators, re- 
spectively. The coupling coefficients are normalized such 
that J2i 9i = l an d individual nuclear spin operators uf 
are assumed to be spin 1/2 for simplicity; g is related to 
the total HF coupling strength Ahf via g = A^p/ J2i 9i- 
We neglect the typically very small nuclear Zeeman and 
nuclear dipole-dipole terma^. For simplicity, we also re- 
strict our analysis to one nuclear species only. These sim- 



plifications will be addressed in more detail in Sec. VI 



The effect of the HF interaction with the nuclear spin 
ensemble is two-fold: The first part of the above Hamil- 
tonian Hff = | (A + S~ + A~ S + ) is a Jaynes-Cummings- 
type interaction which exchanges excitations between the 
QD electron and the nuclei. The second term Hon — 
gA z S z constitutes a quantum magnetic field, the Over- 
hauser field, for the electron spin generated by the nuclei. 
If the Overhauser field is not negligible compared to the 
external Zeeman splitting, it can have a marked effect on 
the current by (de)tuning the hyperfine flip-flops. 



IV. GENERALIZED QUANTUM MASTER 
EQUATION 

Electron transport through a QD can be viewed as 
a tool to reveal the QD's nonequilibrium properties in 
terms of the current-voltage I/V characteristics. From a 
theoretical perspective, a great variety of methods such 
as the scattering matr ix formalism^ and non-equilibrium 
Green's funct ion have been used to explore the I/V 
characteristics of quantum systems that are attached to 
two metal leads. Our analysis is built upon the master 
equation formalism, a tool widely used in quantum op- 
tics for studying the irreversible dynamics of quantum 
systems coupled to a macroscopic environment. 

In what follows, we will employ a projection opera- 
tor based technique to derive an effective master equa- 
tion for the QD system - comprising the QD electron 
spin as well as the nuclear spins - which experiences dis- 
sipation via the electron's coupling to the leads. This 
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dissipation is shown to dynamically depend on the state 
of the nuclear system potentially resulting in feedback 
mechanisms between the two subsystems. We will derive 
conditions which allow for a Markovian treatment of the 
problem and list the assumptions our master equation 
based framework is based on. 



A. Superoperator Formalism - Nakajima-Zwanzig 
Equation 

The state of the global system that comprises the QD 
as well as the environment is represented by the full den- 
sity matrix p(t). However, the actual states of interest 
are the states of the QD which are described by the re- 
duced density matrix ps — Tre [p] , where Tre . . . aver- 
ages over the unobserved degrees of freedom of the Fermi 
leads. We will derive a master equation that governs the 
dynamics of the reduced density matrix p§ using the su- 
peroperator formalism. We start out from the von Neu- 
mann equation for the full density matrix 



p = -i[H (t),p] 



(7) 



where H (t) can be decomposed into the following form 
which turns out to be convenient later on 



H (t) = H (i) + H x (t) + H T . 



(8) 



Here, H (t) = Hz + Hb + g (A z ) t S z comprises the Zee- 
man splitting caused by the external magnetic field via 
Hz and the Hamiltonian of the non-interacting electrons 
in the leads Hb', moreover, the time-dependent expec- 
tation value of the Overhauser field has been absorbed 
into the definition of Hq (t). The HF interaction between 
the QD electron and the ensemble of nuclear spins has 
been split up into the flip-flop term Hg and the Over- 
hauser field Hon, that is Huf = Hon + Hg. The 
term Hi (t) = Haoh (t) + Hg comprises the Jaynes- 
Cummings-type dynamics Hg and fluctuations due to 
deviations of the Overhauser field from its expectation 
value, i.e., Haoh (t) = gSA z S z , where 6A Z = A z - (A z ) t . 

The introduction of superoperators - operators act- 
ing on the space of linear operators on the Hilbert space 
- allows for a compact notation. The von Neumann 
equation is written as p = —iC (t) p, where C (t) = 
£o (t) +C\ (t)+CT is the Liouville superoperator defined 
via C a - = [H a , •]. Next, we define the superoperator V 
as a projector onto the relevant subspace 



Vp (t) = Tr B [p (*)] ® p° B = Ps (t) ® P° B , 



(9) 



where p B describes separate thermal equilibria of the two 
leads whose chemical potentials are different due to the 
bias voltage eV = p^ — pR- Essentially, V maps a density 
operator onto one of product form with the environment 
in equilibrium but still retains the relevant information 
on the system state. The complement of V is Q = 1—V. 

By inserting V and Q in front of both sides of the von 
Neumann equation one can derive a closed equation for 



the projection Vp(t), which for factorized initial condi- 
tion, where Qp (0) = 0, can be rewritten in the form of 
the generalized Nakajima-Zwanzig master equation 

d 

—Tp = -iPCVp 
dt 

rt 

- / dt'VCQfe- l ^ dTQC{T) QCVp(t')ilQ) 
Jo 

which is non-local in time and contains all orders of the 
system-leads coupling^. Here, T denotes the chronolog- 
ical time-ordering operator. Since V and Q are projec- 
tors onto orthogonal subspaces that are only connected 
by £t, this simplifies to 



d_ 

dt 



Vp = -iVCTp- [ dt'V£ T f e - l ti dTQC ^C T V P (*') . 
Jo 

(11) 

Starting out from this exact integro-differential equation, 
we introduce some approximations: In the weak cou- 
pling limit we neglect all powers of Ct higher than two 
(Born approximation) . Consequently, we replace C (t) 
by C (t) — Ct in the exponential of Eqn.( 11 ). Moreover, 



we make use of the fact that the nuclear spins evolve 
on a time-scale that is very slow compared to all elec- 
tronic processes: In other words, the Overhauser field is 
quasi-sta tic on the timescale of single electronic tunnel- 
ing event ^221221. That is, we replace (A Z ) T by (A z ) t in the 
exponential of Eqn.(ll) which removes the explicit time 
dependence in the kernel. By taking the trace over the 
reservoir and using Tre [Vp (<)] = ps (t), we get 



(12) 



PS {t) = -l{C Z +CliF)ps (t) 
rt , 

dTTr B (c T e- l[Ca[t)+Cl{t)]T C T Vp{t-T) 
o ^ 

Here, we also used the relations VCtV = and CbV = 
and switched the integration variable to r = t — t! ' . 
Note that, for notational convenience, we will suppress 
the explicit time-dependence of £o(i) (t) in the following. 
In the next step, we iterate the Schwingcr-Dyson identity 



e ~i(C a +Ci)T _ e -iC T 



(13) 



-l I dr' e -^o{r-r') Lie -i{C +Uy 

'o 



In what follows, we will keep only the first term of this 
infinite series (note that the next two leading terms are 
explicitly calculated in Appendix |A|. In quantum optics, 
this simplification is well known as approximation of in- 
dependent rates of variatiorP^. In our setting it is valid, 
if C\ (t) is small compared to C$ (t) and if the bath cor- 
relation time t c is short compared to the HF dynamics, 
Ahf <C 1/t c . Pictorially, this means that during the cor- 
relation time t c of a tunneling event, there is not sufficient 
time for the Rabi oscillation with frequency g < ^4hf to 
occur. For typical materials^, the relaxation time r c is in 
the range of ~ 10 -15 s corresponding to a relaxation rate 
T c = r^ 1 w 10 s ^eV. Indeed, this is much faster than 
all other relevant processes. In this limit, the equation 
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of motion for the reduced density matrix of the system 
simplifies to 



Ps (t) = -i (Cz + £hf) PS (*) 



(14) 



dr Trt 



(c T e- iC °^jC TPs (t-r)(g>p B ) 



Note, however, that this master equation is not Marko- 



vian as the rate of change of ps (t) still depends on its 
past. Conditions which allow for a Markovian treatment 
of the problem will be addressed in the following. 



B. Markov Approximation 

Using the general relation er iC ° T = e - tHaT Oe lH ° T 
for any operator O, we rewrite Eqn.(14) as 



ps (t) = -i[H z + H HF ,p s (t)} 



H T ( T ),e- lll ° T p s (t-r) 



P°B 



(15) 



In accordance with the previous approximations, we will 
replace e~ lH ° T ps (t — t) e • ° r by ps (t) which is approx- 
imately the same since any correction to H p would be 
of higher order in perturbation theory^lMI. In other 
words, the evolution of ps (t — r) is approximated by 
its unperturbed evolution which is legitimate provided 
that the relevant timescale for this evolution r c is very 
short (Markov approximation). This step is motivated 
by the typically rapid decay of the lead correlations 
functional; the precise validity of this approximation is 
elaborated below. In particular, this simplification disre- 
gards dissipative effects induced by Ht which is valid self- 
consistently provided that the tunneling rates are small 
compared to the dynamics generated by H . 

Moreover, in Eqn.(15) we introduced the tunneling 



Hamiltonian in the interaction picture as Ht (t) 



-iHnr 



» T H T e 



Hot 



For simplicity, we will only consider one 



lead for now and add the terms referring to the second 
lead later on. Therefore, we can disregard an additional 
index specifying the left or right reservoir and write ex- 
plicitly 



H T (r) = J2 r^e-^'M-^ ' d\c ka + h.c. (16) 

k j cr 



Here, the resonances e a (t) are explicitly time-dependent 
as they dynamically depend on the polarization of the 
nuclear spins 



(t) 



± g - (A z 
2 V 



The quantity 



oj = e t (i)-ej. (t) =u Q +g(A z ) t 



(17) 



(18) 



can be interpreted as an effective Zeeman splitting which 
incorporates the external magnetic field as well as the 
mean magnetic field generated by the nuclei. 

Since the leads are assumed to be at equilibrium, their 
correlation functions are given by 



Tr B 4a (r) c fe'*'Ps 



fk 



(19) 



Tr E 



= 5 a , a ,6 k , k ,e ie * T (l- f k ), (20) 



where the Fermi function f k = (1 + exp [j3 (e k — fi)]) 1 
with inverse temperature /3 = l/fcgT gives the thermal 
occupation number of the respective lead in equilibrium. 
Note that all terms comprising two lead creation c£ ct or 
annihilation operators c ka vanish since p a B contains states 
with definite electron number onljE^I. The correlation 
functions are diagonal in spin space and the tunneling 
Hamiltonian preserves the spin projection; therefore only 
co-rotating terms prevail. If we evaluate all dissipative 



terms appearing in Eqn.(15), due to the conservation of 



momentum and spin in Eqn.(19l and Eqn.(20l, only a 
single sum over k, a survives. Here, we single out one 
term explicitly, but all other terms follow analogously. 
We obtain 

Ps (*) = •• -+V /* drC a (r) 4e- iH ° T Ps (* - r) e iH ° T d„, 

( 21 ) 

where the correlation time of the bath t c is determined 
by the decay of the noise correlations 

C a (r)=J2\Tk,a\ 2 fke iMt) - £k)r 

deJ <T (e)e i(e " (t) - e)r . (22) 

Here, we made use of the fact that the leads are macro- 
scopic and therefore exhibit a continuous density of states 
per spin n(e). On top of that, we have introduced the 
spectral density of the bath as 



J a (e) = D a (e) / (e) , 



(23) 



where D a (e) = n(e) \T„ (e)| is the effective density of 
states. The Markovian treatment manifests itself in a 
self-consistency argument: We assume that the spec- 
tral density of the bath J a (e) is flat around the (time- 
dependent) resonance e a (t) over a range set by the char- 
acteristic width Ttj. Typically, both the tunneling ma- 
trix elements T a (e) as well as the density of states n (e) 
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FIG. 3. (color online). Fermi function for finite temperature 
(dashed blue line) and in the limit T = (solid blue line). 
The absolute value of the derivative of the Fermi function 
/'(e) (dotted orange line for finite temperature) is maximized 
at the chemical potential fj, and tends to a delta function in 
the limit T — > 0. The Markovian description is valid provided 
that the Fermi function is approximately constant around 
the resonances e CT (t) on a scale of the width of these reso- 
nances, schematically shown in red (solid line for e a (t) < fi 
and dashed line for e CT (t) > fi). 



are slowly varying functions of energy. In the so-called 
wide-band limit the effective density of states D a (e) is 
assumed to be constant so that the self-consistency argu- 
ment will exclusively concern the behaviour of the Fermi 
function / (e) which is intimately related to the temper- 
ature of the bath T. Under the condition, that J a (e) 
behaves fiat on the scale T^, it can be replaced by its 
value at e a (t), and the noise correlation simplifies to 



C CT (r) = J a {e a (i))e 



i<L a (t)r 



de t 



Using the relation 



irS (r) — il 



(24) 



(25) 



with P denoting Cauchy's principal value, we find that 
the Markov approximation Re [C a (r)] oc S (r) is ful- 
filled provided that the self-consistency argument holds. 
This corresponds to the white-noise limit where the 
correlation-time of the bath is r c = 0. Pictorially, 
the reservoir has no memory and instantaneously re- 
laxes to equilibrium. We can then indeed replace 



-iHnr 



PS (t 



t e 



by ps {t) and extend the integra- 



tion in Eqn.(15l to infinity, with neglibile contributions 
due to the rapid decay of the memory kernel. In the fol- 
lowing, we will derive an explicit condition for the self- 
consistency argument to be satisfied. 

Let us first consider the limit T — 0: As schematically 
depicted in Fig. [3j in this case / (e) behaves perfectly 
flat except for e — fj, where the self-consisteny argument 



is violated. Therefore, the Markovian approximation is 
valid at T = given that the condition \e a (t) — fi\ 3> Td 
is fulfilled. In this limit, all tunneling rates are constant 
over time and effectively decoupled from the nuclear dy- 
namics. Note that for the observation of electronic su- 
perradiance it will be sufficient to restrict oneself to this 
case. 

For a more general analysis, we now turn to the case of 
finite temperature T > 0. We require the absolute value 
of the relative change of the Fermi function around the 
resonance e a (t) over a range of the characteristic width 
Td to be much less than unity, that is 



9f(e) 




de 





r d « i. 



(26) 



An upper bound for the first factor can easily be obtained 
as this quantity is maximized at the chemical potential [i, 
for all temperatures. Evaluating the derivative at e a (i) = 
/j, results in the compact condition 



r d « 4k B T. 



(27) 



Thus, finite temperature T > washes out the rapid 
character of / (e) at the chemical potential y, and, pro- 
vided that Eqn.([27]) is fulfilled, allows for an Markovian 
treatment. 

Two distinct mechanisms contribute to the width T^: 
dissipation due to coupling to the leads and the effect 
of Hi (t), because both of them have been neglected 
sclf-consistently in the memory kernel when going from 
Eqn.(|TT|) to Eqn.(14|. Typically, the tunneling rates are 
of the order of ~ 5 — 20/ieV, depending on the trans- 
parency of the tunnel-barrier. Regarding the contribu- 
tion due to Hi (t), we first consider two limits of partic- 
ular importance: For a completely mixed state the fluc- 
tuation of the nuclear field around its zero expectation 
value is of the order of ~ Anp/y/N s» 0.1 fieV. In con- 
trast, for a fully polarized state these fluctuations can be 
neglected whereas the effective strength of the flip-flop 
dynamics is ~ Ahf / V~N as well. Therefore, in both lim- 
its considered here, the dominant contribution to Td is 
due to the coupling to the leads and the self-consistency 
condition could still be met with cryostatic temperatures 
ksT > 10 /xeV, well below the orbital level spacing. How- 
ever, we note that in the course of a superradiant evo- 
lution, where strong correlations among the nuclei build 
up, the dominant contribution to Td may come from the 
flip-flop dynamics, which are Ahf/4 ~ 25 /ieV at max- 
imum for homogeneous coupling. For realisitic condi- 
tions, though, this effect is significantly reduced, as will 
be demonstrated in our simulations in Sec. IVII 



C. General Master Equation for Nuclear Spin 
Assisted Transport 

Assuming that the self-consistency argument for a 
Markovian treatment is satisfied, we now apply the fol- 
lowing modifications to Eqn.(15): First, we neglect level 



shifts due to the coupling to the continuum states which 
can be incorporated by replacing the bare frequencies 
e a {t) with renormalized frequencies. Second, one adds 
the second electron reservoir that has been omitted in 
the derivation above. Lastly, one performs a suitable 
transformation into a frame rotating at the frequency 
e = (e^ + ei) /2 leaving all terms invariant but changing 

H z from H z = efdldf + £4. 4^4. to H z = luqS z . Af- 
ter these manipulations one arrives at the central master 
equation as stated in Eqn.([2| where the tunneling rates 

With a a (t) = E,=L,fl«^ (*), Pa (t) = E X =L,R& X) (*) 

and 



og (t) 



n x (e CT (t)) (e ff (t)) [1 - / x (e CT (i))] 



^-Ai = % (e CT (t)) T<*> (e CT (<)) / a (e ff (*)) 

govern the dissipative processes in which the QD system 
exchanges single electrons with the leads. The tunneling 
rates, as presented here, are wid ely used in nanostruc- 
ture quantum transport problema^ESEl However, in 
our setting they are evaluated at the resonances e a (t) 
which dynamically depend on the polarization of the nu- 
clear spins; see Eqn.(17). Note that Eqn.Q incorporates 
finite temperature effects via the Fermi functions of the 
leads. This potentially gives rise to feedback mechansims 
between the electronic and the nuclear dynamics, since 
the purely electronic diffusion markedly depends on the 
nuclear dynamics. 

Since Eqn.Q marks our first main result, at this point 
we quickly reiterate the assumptions our master equation 
treatment is based on: 

• The system-lead coupling is assumed to be weak 
and therefore treated perturbatively up to second 
order (Born-approximation). 

• In particular, the tunneling rates are small com- 
pared to the effective Zeeman splitting u>. 

• Level shifts arising from the coupling to the con- 
tinuum states in the leads are merely incorporated 
into a redefinition of the QD energy levels t a (t). 

• There is a separation of timescales between 
electron-spin dynamics and nuclear-spin dynamics. 
In particular, the Overhauser field g(A z ) t evolves 
on a timescale that is slow compared to single elec- 
tron tunneling events. 

• We have applied the approximation of independent 
rates of variation: If the HF dynamics generated 
by Hi (t) = H s + i?AOH (t) is (i) sufficiently weak 
compared to Hq and (ii) slow compared to the cor- 
relation time of the bath t c , their effect in the mem- 
ory kernel of the master equation can be neglected. 
The latter holds for ^4hf t c "C 1. Note that the flip- 
flop dynamics can become very fast as correlations 



among the nuclei build up culminating in a maxi- 
mum coupling strength of Ahf /4 for homogeneous 
coupling. This potentially drives the system into 
the strong coupling regime where condition (i), that 
is u) 3> H-Hi it) ||, might be violated. However, un- 
der realistic conditions of inhomogeneous coupling 
this effect is significantly reduced. 

• The effective density of states D a (e) = 
n (e) \T a (e)| 2 is weakly energy-dependent (wide- 
band limit). In particular, it is flat on a scale of 
the characteristic widths of the resonances. 

• The Markovian description is valid provided that 
either the resonances are far away from the chem- 
ical potentials of the leads on a scale set by the 
characteristic widths of the resonances or the tem- 
perature is sufficiently high to smooth out the rapid 
character of the Fermi functions of the leads. This 



condition is quantified in Eqn.(27) 



In summary, we have derived a Quantum master equa- 
tion describing electronic transport through a single QD 
which is collectively enhanced due to the interaction with 
a large ancilla system, namely the nuclear spin ensemble 
in the host environment. Eqn.([2| incorporates two major 
intriguing features both of theoretical and experimental 
relevance: Due to a separation of timescales, only the 
electronic subsystem experiences dissipation with rates 
that depend dynamically on the state of the ancilla sys- 
tem. This non-linearity gives rise to feedback mecha- 
nisms between the two subsystems as well as hysteretic 
behavior. Moreover, the collective nature of the HF inter- 
action offers the possibility to observe intriguing coherent 
many-body effects. Here, one particular outcome is the 
occurence of electronic superradiance, as will be shown 
in the remainder of this paper. 

Note that in the absence of HF interaction between the 
QD electron and the proximal nuclear spins, i.e., in the 
limit g — ► 0, our results agree with previous theoretical 
studies^. 



V. ELECTRONIC SUPERRADIANCE 

Proceeding from our general theory derived above, this 
section is devoted to the prediction and analysis of super- 
radiant behavior of electrons tunneling through a single 
QD in the Coulomb-blockade regime; see Fig. [T] for the 
scheme of the setup. 

We note that, in principle, an enhancement seen in 
the leakage current could simply arise from the Over- 
hauser field dynamically tuning the hyperfine flip-flops. 
However, we can still ensure that the measured change 
in the leakage current through the QD is due to coop- 
erative emission only by dynamically compensating the 
Overhauser field. This can be achieved by applying a 
time dependent magnetic or spin-dependent AC Stark 
field such that H comp (t) 



-g{A z ) t S z which will be 
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done in most of our simulations below to clearly prove 
the existence of superradiant behaviour in this setting. 
Consequently, in our previous analysis H (t) is replaced 
by H = Ho (i) - g (A z ) t S z = H z + H B so that the po- 
larization dependence of the tunneling rates is removed 
and we can drop the explicit time-dependence of the res- 
onances e a (t) — > e a . Under this condition, the master 
equation for the reduced system density operator can be 
written as 

p S (t) = ~i [UJ S Z + tf HF + tfcomp (t) , PS (t)] (28) 



■ E & 



daPS (t) 4 - 2 { d l d a,PS (t)} 



dtps (*) d a - - {d a dl,p s (t)} 



In accordance with our previous considerations, in this 
specific setting the Markovian treatment is valid provided 
that the spectral density of the reservoirs varies smoothly 
around the (time-independent) resonances e a on a scale 
set by the natural widths of the level and the fluctuations 
of the dynamically compensated Overhauser field. More 
specifically, throughout the whole evolution the levels are 
assumed to be far away from the chemical potentials of 
the reservoir d 39 ! 40 !, for an illustration see Fig. [3] In this 
wide band limit, the tunneling rates a a , /3 CT are indepen- 
dent of the state of the nuclear spins. The master equa- 
tion is of Lindblad form which guarantees to preserve the 
complete positivity and the hermiticity of the density ma- 
trix. Eqn.(28) agrees with previous theoretical results^ 



except for the appearance of the collective HF interaction 
between the QD electron and the ancilla system in the 



Hamiltonian dynamics of Eqn.(28). 

To some extent, Eqn.(p8| bears some similarity with 



the quantum theory of the laser. While in the latter the 
atoms interact with bosonic reservoirs, in our transport 

- b w 

and emits fermionic particlej 45 * 46 l 

If the HF dynamics are the slowest timescale in the 



setting the QD is pumped b y the nuclear spin ensemble 



problem, Eqn.(28 1 can be recast into a form which makes 
its superradiant character more apparent. In this case, 
the system is subject to the slaving principled The dy- 



namics of the whole system follow that of the subsystem 
with the slowest time constant allowing to adiabatically 
eliminate the electronic QD coordinates and to obtain 
an effective equation of motion for the nuclear spins. In 
this limit, the Overhauser field is much smaller than the 
Zeeman splitting so that a dynamic compensation of the 
OH can be disregarded for the moment. For simplic- 
ity we consider a transport setting in which only four 
tunneling rates are different from zero, see Fig. [1] The 
QD can be recharged from the left and the right lead, 
but only electrons with spin projection a =f can tunnel 
out of the QD into the right lead. We define the total 
recharging rate j3 



pi L ^ and 



for notational convenience unambiguously set a = a 



(-R) 



First, we project Eqn.(28) onto the populations of the 




FIG. 4. (color online). The electronic QD system in the 
local moment regime after the adiabatic elimination of the 
|0) level including the relevant dissipative processes. Within 
the effective system (box) we encounter an effective decay 
term and an effective pure dephasing term, given by the rates 
7 and F respectively. This simplification is possible for fast 
recharging of the QD, i.e., /3 2> 01. 



electronic levels and the coherences in spin space accord- 
ing to p mn = (m\ ps \n), where m, n — 0, f, |. This yields 



Poo = "Pit ~ ^<°oo ( 29 ) 

Pn = -*2 t^ 2 ' Prd ~ *| ^ A ~ p ^ ~ p^ a+ ) ~ a ^tt + /Voo (30) 

PU = +*§ W,Pu\ ~ *f ( A+ Ptl ~ PlfA~) + ft/loo (31) 

pu = - iuj oPn - *2 ^ A *^ + Pn AZ ) ~ *| i A ~Pu - Pn A ~) - ~^Pn- (32) 



We can retrieve an effective master equation for the ciently strong recharging rate, that is in the limit f3 3> a, 
regime in which on relevant timescales the QD is al- which can be implemented experimentally by making the 
ways populated by an electron. This holds for a suffi- left tunnel barrier more transparent than the right one. 
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Then, the state |0) is populated negligibly throughout 
the dynamics and can be eliminated adiabatically accord- 
ing to poo » % Pff. In analogy to the Anderson impurity 
model, in the following this limit will be referred to as lo- 
cal moment regime. The resulting effective master equa- 
tion reads 



ps = —i [uoS z + -Hhf i Ps] 



(33) 



+7 

+r 



S-psS + --{S + S-,p s } 



S Z PS S Z - ~ps 



where 7 = ^ a is an effective decay rate and r = ^aan 
effective dephasing rate. This situation is schematized 
in Fig. [4] The effective decay (dephasing) describes pro- 
cesses in which the QD is recharged with a spin down 
(up) electron after a spin up electron has tunneled out of 
the QD. 

In the next step we aim for an effective description 
that contains only the nuclear spins: Starting from a fully 
polarized state, SR is due to the increase in the operative 
HF matrix element (A + A~) . The scale of the coupling 
is set by the total HF coupling constant A RF = g V . g t . 
For a sufficiently small relative coupling strengtfP^ 



where 



e = A H f/(2A), 



A = \a/2 + iuj \ 



(34) 



(35) 



the electron is predominantly in its \l) spin state and 
we can project Eqn.(33) to the respective subspace. As 
shown in detail in Appendix [B] in this limit the re- 
duced master equation for the nuclear density operator 
p = Tr e | [ps] is given by Eqn. ([3]), where the effective co- 
efficients read 



g 2 g 

4A 2 ' 
4A 2 ' 



(36) 
(37) 



This master equation is our second main result. In an 
optical setting, it has previously been predicted theoreti- 
cally to exhibit strong SR signatures^. Conceptually, its 
superradiant character can be understood immediately 
in the ideal case of homogeneous coupling in which the 
collective state of all nuclear spins can be described in 
terms of Dicke states \J,m): The enhancement of the 
HF interaction is directly associated with the transition 
through nuclear Dicke states | J, m) , m <C J- In this ideal- 
ized setting, the angular momentum operator I = \/N A 
of the nuclear spin ensemble obeys the SU(2) Lie al- 
gebra, from which one can deduce the ladder operator 
relation I~ \J, to) = */ J (J + 1) — Tn (to — 1) \J,m— 1). 
This means that, starting from an initially fully polarized 



state I J = A/2, m = A/2) , the system cascades down the 
Dicke-ladder with an effective rate 



r, 



N 



(N/2 + to) (A/2 - to + 1) 



(38) 



since, according to the first term in Eqn.([3|, the popula- 
tions of the Dicke states evolve as 



-^(A/2 + to) ( A/2 - m + 1) Mr, 



(39) 



+ ^ (A/2 + TO + 1) (A/2 - TO) p m +l,m+l. 



While the effective rate is T 



N/2-+N/2-1 = c r at the very 
top of the the ladder it increases up to ^\ m \<^N/2 ~ i[N 
at the center of the Dicke ladder. This implies the char- 
acteristic intensity peaking as compared to the limit of 
independent classical emitters the emission rate of which 
would be T cl = f A t = % (N/2 + to). 

However, there is also a major difference compared to 
the superradiant emission of photons from atomic en- 
sembles: In contrast to its atomic cousin, the prefactor 
c r /N oc 1/A 2 is A-dependent, resulting in an overall 
time of the SR evolution (to) which increases with N. By 
linearizing Eqn. ( 38 ) for the beginning of the superradiant 

c r (s - 



evolution^ as r m _ > . m _ 1 c r (s + 1), where s = N/2 — to 
gives the number of nuclear flips, one finds that the first 
flip takes place in an average time c" 1 , the second one 
in a time (2c r )~ 1 and so on. The summation of all these 
elementary time intervals gives an upper bound estimate 
for the process duration till the SR peaking as 



< 



1 



1 



N/2 



21n(A/2) 



(40) 



which, indeed, increases with the number of emitters as 
~ Aln(A), whereas one obtains (tn) ~ \n(N)/N for or- 
dinary superradiance^l Accordingly, in our solid-state 
system the characteristic SR peak appears at later times 
for higher N. The underlying reason for this difference is 
that in the atomic setting each new emitter adds to the 
overall coupling strength, whereas in the central spin set- 
ting a fixed overall coupling strength Ahf is distributed 
over an increasing number of particles. Note that in 
an actual experimental setting N is not a parameter, of 
course. For our theoretical discussion, though, it is con- 
venient to fix the total HF coupling strength ^4hf and to 
extrapolate from our findings to an experimentally rele- 
vant number of nulear spins A. 

For large relative coupling strength e » 1 the QD 
electron saturates and superradiant emission is capped 
by the decay rate a/2, prohibiting the observation of a 
strong intensity peak. In order to circumvent this bottle- 
neck regime, one has to choose a detuning uio such that 
< e < 1. However, to realize the spin-blockade regime, 
where the upper spin manifold is energetically well sepa- 
rated from the lower spin manifold, the Zeeman splitting 
has to be of the order of u>o ~ ^4hf which guarantees 
e < 1. In this parameter range, the early stage of the 
evolution - in which the correlation buildup necessary 
for SR takes place^- is well described by Eqn. ([3]). 
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In reality, the inhomogeneous nature ^ const) of 
the collective operators leads to dephasing between 
the nuclei, possibly preventing the phase d emission nec- 
essary for the observation of SR 24 25 36 37 . The inhomoge- 
neous part of last term in Eqn. ( 3 ) - the electron's Knight 
field - causes dephasing^! oc gryVar (<7i)/2, possibly lead- 
ing to symmetry reducing transitions J — > J — 1. Still, 
it has been shown that SR is also present in realistic 
inhomogeneous systems^, since the system evolves in a 
many-body protected manifold (MPM) : The second term 
m energetically separates different total nuclear 

spin-J manifolds protecting the correlation build-up for 
large enough e. 

The superradiant character of Eqn. Q suggests the ob- 
servation of its prominent intensity peak in the leakage 
current through the QD in the spin-blockade regime. We 
have empl oyed the method of Full-Counting-Statistics 
(FCSpi^S! j n orc j er to obtain an expression for the cur- 
rent and find (setting the electron's charge e = 1) 



I (t) = ap n 



Pi Poo- 



(41) 



This result is in agreement with previous theoretical find- 
ings: The current through the device is completely deter- 
mined b y the oc cupation of the levels adjacent to one of 
the lead d 39 * 42 ^ 9 -! The first term describes the accumula- 
tion of electrons with spin a =t in the right lead, whereas 
the second term describes electrons with a =4- tunneling 
from the right lead into the QD. As done before^, we 
take the ratio of the maximum current to the initial cur- 
rent (the maximum for independent emitters) I COO p/Iind 
as our figure of merit: a relative intensity peak height 
-^coopMnd > 1 indicates cooperative effects. One of the 
characteristic features of SR is that this quantity scales 
linearly with the number of spins N. 



In the local-moment regime, described by Eqn. (33 1, 
the expression for the current simplifies to I (t) — (T— 

P ! f ) /P)a (S+S-) t cx (S+S-) t showing that it is directly 
proportional to the electron inversion. This, in turn, in- 
creases as the nuclear system pumps excitations into the 
electronic system. A compact expression for the relation 
between the current and the dynamics of the nuclear sys- 
tem can be obtained immediately in the case of homoge- 
neous coupling 



d (s + s-) t = -f t (i*) t -i(s + s-) t 



(It 



(42) 



Since the nuclear dynamics are in general much slower 
than the electron's dynamics, the approximate solution 
of this equation is (S + S~) t ~ — (I z ) t /j- As a con- 
sequence, the current / (t) is proportional to the time- 
derivative of the nuclear polarization 



(43) 



Still, no matter how strong the cooperative effects are, 
on a timescale of single electron tunneling events, the 
electrons will always be emitted antibunched, since in the 



strong Coulomb-blockade regime the QD acts as a single- 
electron emitter^. Typically, the rate for single-electron 
emission events is even below the single tunneling rate a 
due to the spin-blockade. On electronic timescales ~ 1/a, 
the SR mechanism manifests in lifting this blockade; as 
argued above, the efficiency of this process is significantly 
enhanced by collective effects. 

Before we proceed with an in-depth analysis of the 
current I(t), we note that an intriguing extension 
of the present work would be the study of fluctua- 
tions thereof. Insights into the statistics of the cur- 
rent could be obtained by analyzing two-time corre- 
lation functions such as (rif(t + T)nf(t)), where — 
d^df. This can conveniently be done via the Quan- 
tum Regression Theorem^ which yields the formal re- 
sult (n t (t + T)n t (t)) = Tr s [n t e Wr (n t p s {t))] . Here, W 
denotes the Liouvillian governing the system's dynamics 
according to p s = Wps and Trg [. . . ] refers to the trace 
over the system's degree of freedoms. This procedure 
can be generalized to higher order correlation functions 
and full evaluation of the current statistics might reveal 
potential connections between current fluctuations and 
cooperative nuclear dynamics. 



VI. ANALYSIS AND NUMERICAL RESULTS 

A. Experimental Realization 

The proposed setup described here may be realized 
with state-of-the-art experimental techniques. First, the 
Markovian regime, valid for sufficiently large bias eV , is 
realized if the Fermi functions of the leads are smooth on 
a scale set by the natural widths of the levels and residual 
fluctuations due to the dynamically compensated Over- 
hauser field. Since for typical materials the hyperfine 
coupling constant is Ahf = 1 — 100/ieV^ and tunneling 
rates are typically of the order of ~ 10 pe\ffi, this does 
not put a severe restriction on the bias voltage whic h are 
routinely in the range of hundreds of fiV or m V 18 * 19 [ Sec- 
ond, in order to tune the system into the spin-blockade 
regime, a sufficiently large external magnetic field has to 
be applied. More precisely, the corresponding Zeeman 
splitting wo energetically separates the upper and lower 
manifolds in such a way that the Fermi function of the 
right lead drops from one at the lower manifold to zero 
at the upper manifold. Temperature smeares out the 
Fermi function around the chemical potential by approx- 
imately ~ kgT. Accordingly, with cryostatic tempera- 
tures of ksT ~ 10 peV being routinely realized in the 
labS^l, this condition can be met by applying an exter- 
nal magnetic field of ~ 5 — 10T which is equivalent to 
lu « 100 - 200^eV in GaAfP^U. Lastly, the charging en- 
ergy U, typically ~ 1 — 4meV |9J1 ^, sets the largest energy 
scale in the problem justifying the Coulomb-blockade 
regime with negligible double occupancy of the QD pro- 
vided that the chemical potential of the left lead is well 
below the doubly occupied level. Lastly, we note that 
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FIG. 5. (color online). Typical time-evolution of the nor- 
malized current for homogeneous coupling under dynamical 
compensation of the Overhauser field and a relative coupling 
strength of e = 0.5, shown here for TV = 60 and TV = 100 
nuclear spins. The characteristic feature of superradiance, a 
pronounced peak in the leakage current proportional to TV, is 
clearly observed. 



similar setups to the one proposed here have prev iously 
been realized experimentally by Hanson et al! 34 * 35 *. 

Proceeding from these considerations, we will now 
show by numerical simulation that an SR peaking of sev- 
eral orders of magnitude can be observed for experimen- 
tally relevant parameters in the leakage current through 
a quantum dot in the spin-blockade regime. We will first 
consider the idealized case of homogeneous coupling for 
which an exact numerical treatment is feasible even for a 
larger number of coupled nuclei. Then, we will continue 
with the more realistic case of inhomogeneous coupling 
for which an approximative scheme is applied. Here, we 
will also study scenarios in which the nuclear spins are 
not fully polarized initially. Finally, we will discuss in- 
trinsic nuclear dephasing effects and undesired cotunnel- 
ing processes which have been omitted in our simulations. 



B. Superradiant Emission of Electrons 

1. Idealized Setting 

The homogeneous case allows for an exact treatment 
even for a relatively large number of nuclei as the system 
evolves within the totally symmetric low-dimensional 
subspace {| J, m) ,m = — J, . . . , J}. Starting from a fully- 
polarized state, a strong intensity enhancement is ob- 
served; typical results obtained from numerical simula- 
tions of Eqn.(28) are depicted in Fig. [5] for N = 60 and 
N = 100 nuclear spins. The corresponding relative peak 
heights display a linear dependence with N, cf. Fig. |6j 
which we identify as the characteristic feature of super- 
radiance. Here, we have used the numerical parameters 



1 HF 



1, uj — 1 and a 



L) 



0.1 in 



units of ~ 100 /ieV, corresponding to a relative coupling 
strength e = 0.5. 

Before we proceed, some further remarks on the dy- 
namic compensation of the Overhauser field seem ap- 
propriate: We have merely introduced it in our analysis 
in order to provide a clear criterion for the presence of 
purely collective effects, given by / C oop/^ind > 1. In other 
words, dynamic compensation of the Overhauser field is 
not a necessary requirement for the observation of collec- 
tive effects, but it is rather an adequate tool to display 
them clearly. From an experimental point of view, the 
dynamic compensation of the Overhauser field might be 
challenging as it requires accurate knowledge about the 
evolution of the nuclear spins. Therefore, we also present 
results for the case in which the external magnetic field is 
constant and no compensation is applied. This markedly 
relaxes the experimental challenges. Here, we can distin- 
guish two cases: Depending on the sign of the HF cou- 
pling constant ^4hf, the time-dependence of the effective 
Zeeman-splitting uj can either give rise to an additional 
enhancement of the leakage current (Ahf > 0) or it can 
counteract the collective effects (^4hf < 0). As shown in 
Fig. [6j this sets lower and upper bounds for the observed 
enhancement of the leakage current. 

In Fig. [| we also compare the results obtained for dy- 
namic compensation of the Overhauser field to the ideal- 
ized case of perfect compensation in which the effect of 
the Overhauser field is set to zero, i.e., i?oH = gA z S z = 
0. Both approaches display the same features justifying 
our approximation of neglecting residual (de)tuning ef- 
fects of the dynamically compensated Overhauser field 
w.r.t. the external Zeeman splitting ljq. This will also 
be discussed in greater detail below. 



2. Beyond the Idealized Setting 

We now turn to the more realistic case of inhomoge- 
neous coupling which in principle could prevent the phas- 
ing necessary for SR. However, as shown below, SR is 
still present in realistic inhomogeneous systems. In con- 
trast to the idealized case of homogeneous coupling, the 
dynamics cannot be restricted to a low-dimensional sub- 
space so that an exact numerical treatment is not fea- 
sible due to the large number of nuclei. We therefore 
use an approximate approach which has previously been 
shown to capture the effect of nuclear spin coherences 
whil e allo wing for a numerical treatment of hundreds of 
spin d 22 * 2 ^ . For simplicity, we restrict ourselves to the lo- 
cal moment regime in which the current can be obtained 
directly from the electron inversion / (t) oc (S + S~) t . By 
Eqn.( |33| , this expectation value is related to a hierachy 
of correlation terms involving both the electron and nu- 
clear spins. Based on a Wick type factorization scheme, 
higher order expressions are factorized in terms of the 
covariance matrix 7^ = ((7 2 + o\7) and the "mediated co- 
variance matrix" 7^ = (af S z aj) . For further details, 
see RefsP^^I. Moreover, the coupling constants gj have 
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FIG. 6. (color online). Ratio of the maximum current to 
the initial current I coop /Ii n d as a function of the number of 
nuclear spins N for homogeneous coupling and a relative cou- 
pling strength of e = 0.5: Results for perfect compensation 
(dashed line) are compared to the case of dynamic compensa- 
tion (dotted line) of the Overhauser field (OHC). Simulations 
without compensation of the Overhauser field set bounds for 
the enhancement of the leakage current, depending on the 
sign of the HF coupling constant Ahf ; solid and dash-dotted 
line for ^4hf > and Ahf < 0, respectively. 



been obtained from the assumption of a two-dimensional 
Gaussian spatial electron wavefunction of width y/~N /2. 
Specifically, we will present results for two sets of nu- 
merical parameters, corresponding to a relative coupling 
strength of e = 0.5, where Ahf = 1, = 1, 7 = 0.1 and 
r = 0.08, and e = 0.55 with A HF = 1, w = 0.9, 7 = 0.1 
and T = 0.067. 

As shown in Fig. [7] and [5J the results obtained with 
these methods demonstrate clear SR signatures. In com- 
parison to the ideal case of homogeneous coupling, the 
relative height is reduced, but for a fully polarized ini- 
tial state we still find a linear enhancement I C oop/Imd ~ 
0.043iV (e = 0.5); therefore, as long as this linear depen- 
dence is valid, for typically TV s» 10 5 — 10 6 a strong in- 
tensity enhancement of several orders of magnitude is 
predicted (~ 10 3 - 10 4 ). 

If the initial state is not fully polarized, SR effects are 
reduced: However, when starting from a mixture of sym- 
metric Dicke states | J, J) with polarization p = 80(60)%, 
we find that the linear N dependence is still present: 
/coopMnd ~ 0.0075(0.0028)/V (e = 0.5), i.e., the scal- 
ing is about a factor of ~ 5(15) weaker than for full 
polarizatiorPn Still, provided the linear scaling holds 
up to an experimentally realistic number of nuclei N 
10 5 — 10 6 , this amounts to a relative enhancement of the 
order of / COO p// in d ~ 10 2 - 10 3 . 

In our simulations we have self-consistently verified 
that the fluctuations of the Overhauser field, defined via 
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FIG. 7. (color online). Typical time-evolution of the normal- 
ized current for inhomogeneous coupling, shown here for up 
to N = 13 2 nuclear spins and a relative coupling strength 
e — 0.55. Compared to the idealized case of homogeneous 
coupling, the SR effects are reduced, but still clearly present. 
A Gaussian spatial electron wave function has been assumed 
and the Overhauser field is compensated dynamically. 
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FIG. 8. (color online). Ratio of the maximum current to the 
initial current Icoop/^md as a function of the number of nuclear 
spins TV for relative coupling strengths e = 0.5 and e = 0.55: 
Results for inhomogenous coupling. The linear dependence 
is still present when starting from a nuclear state with finite 
polarization p = 0.8. 



are indeed small compared to the external Zeeman split- 
ting wo throughout the entire evolution. This ensures 
the validity of our perturbative approach and the real- 
ization of the spin-blockade regime. From atomic su- 
perradiance it is known that in the limit of homogeneous 
coupling large fluctuations can build up, since in the mid- 
dle of the emission process the density matrix becomes a 
broad distribution over the Dicke stated^. Accordingly, 
in the idealized, exactly solvable case of homogeneous 
coupling we numerically find rather large fluctuations of 
the Overhauser field; as demonstrated in Fig. |9j this 
holds independently of TV. In particular, for a relative 
coupling strength e = 0.5 the fluctuations culminate in 
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FIG. 9. (color online). Fluctuations of the Overhauser field 
relative to the external Zeeman splitting uiq . In the limit of ho- 
mogeneous HF coupling, strong fluctuations build up towards 
the middle of the emission process (red line, e = 0.5). For 
inhomogeneous coupling this build-up of fluctuations is hin- 
dered by the dephasing between the nuclear spins, resulting 
in considerably smaller fluctuations: The value of the Over- 
hauser fluctuations is shown at the time of the SR peak t max 
for e — 0.5 (orange squares) and e = 0.55 (green diamonds). 
The Overhauser fluctuations reach a maximum value later 
than tmax, see blue dots for e = 0.5. For independent homo- 
geneously coupled nuclear spins, one can estimate the fluctu- 
ations via the Binominal distribution (black line). 



max [Aoh] A^o ~ 0.35. However, in the case of inhomo- 
geneneous HF coupling the Overhauser field fluctuations 
are found to be smaller as the build-up of these fluctua- 
tions is hindered by the Knight term causing dephasing 
among the nuclear spins. As another limiting case, we 
also estimate the fluctuations for completely independent 
homogeneously coupled nuclear spins vi a the Binominal 
distribution as max [Aoh] ~ 0.5Auf/\/N^. Moreover, 
we have also ensured self-consistently the validity of the 
perturbative treatment of the flip-flop dynamics; that is, 
throughout the entire evolution, even for maximum oper- 
ative matrix elements (A + A~) t , the strength of the flip- 
flop dynamics ||-ffff|| was still at least five times smaller 
than wo- 

Initially, the HF mediated superradiance dynamics is 
rather slow, with its characteristic time scale set by c" 1 ; 
for experimentally realistic parameters - in what follows 
we use the parameter set (e — 0.5, awlO /ieV, N w 10 5 ) 
for numerical estimates - this corresponds to c" 1 rs 10 /is. 
Based on fits as shown in Fig. [TUJ we then estimate for 
the SR process duration {to) ~ 50cr* « 500 /is which 
is still smaller than recently reported^] nuclear decoher- 
ence times of ~ 1 ms. Therefore, it should be possible 
to observe the characteristic enhancement of the leakage 
current before the nuclear spins decohere. 

Accordingly, in the initial phasing stage, the HF me- 
diated lifting of the spin-blockade is rather weak result- 
ing in a low leakage current, approximatively given by 
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FIG. 10. (color online). Total time till the observation of 
the characteristic SR peaking i max for e = 0.5 (blue dots) and 
e = 0.55 (orange squares). Based on Eqn.(40 1, logarithmic fits 
are obtained from which we estimate i max for experimentally 
realistic number of nuclear spins N w 10 5 . 



I (t = 0) m e 2 a/N. Therefore, the initial current due 
to HF processes is inversely proportional to the number 
of nuclear spins N. However, as correlations among the 
nuclei build up, the HF mediated lifting becomes more ef- 
ficient culminating in a maximum current of I max ~ e 2 a, 
independent of N. For realistic experimental values - 
also taking into account the effects of inhomogeneous 
HF coupling and finite initial polarization p rj 0.6 - 
we estimate the initial (maximum) leakage current to be 
of the order of I (t = 0) w 6fA (7 max w 10 pA). Leak- 
age currents in this range of magnitudes have already 
been detected in single QD spin-filter experi ments^ as 
well as double QD Pauli-blockade experiment d 15 | 16 | 18 | :L9 !; 
here, leakage currents below 10 fA and 150 fA, respec- 
tively, h ave b een attributed explicitly to other spurious 
processe d 18 * 34 ^. Among others, these will be addressed in 
greater detail in the following. 

In our simulations we have disregarded species inho- 
mogeneities in the nuclear Zeeman energies. In princi- 
ple, these are large enough to cause additional dephas- 
ing between the nuclear spins, similar to the inhomoge- 
neous Knight field. However, this dephasing mechanism 
only applies to nuclei of different Zeeman energies, that 
is nuclei which belong to different species^. This leads 
to two or three mutually decohered subsystems each of 
which is described by our theory. Moreover, we have ne- 
glected the dipolar and quadrupolar interactions among 
the nuclear spins. First, the latter is absent for nuclear 
spin 7 = 1/2 (CdSe QDs) or strain-free QDsP. Sec- 
ond, the nuclear dipole-dipole interaction can cause dif- 
fusion and dephasing processes. Diffusion processes that 
can change A z are strongly detuned and therefore of mi- 
nor importance, as corro borat ed by experimentally mea- 
sured spin diffusion rate d 29 * 31 ^. Resonant processes such 
as oc I?H can lead to dephasing similar to the inho- 
mogenous Knight shift. This competes with the phasing 
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necessary for the observation of SR as expressed by the 
first term in Eqn.Q. The SR process is the weakest at 
the very beginning of the evolution where we estimate 
its strength as c™ ln s» 10 fieV/N. An upper bound for 
the dipole-dipole interaction in GaAs has been given in 
RefpSJas ~ 10 -5 /ieV, in agreement with values given in 
RefsP^H. Therefore, the nuclear dipole-dipole interac- 
tion can safely be neglected for N < 10 5 . In particu- 
lar, its effect should be further reduced for highly polar- 
ized ensembles. Moreover, as argued above due to the 
presence of the MPM-term in Eqn. ^ and demonstrated 
by our simulations, the observation of SR is even robust 
against dephasing caused by the much stronger Knight 
field. 

Our transport setting is tuned into the sequential tun- 
neling regime and therefore we have disregarded cotun- 
neling processes which are fourth order in Ht- In prin- 
ciple, cotunneling processes could lift the spin-blockade 
and add an extra contribution to the leakage current that 
is independent of the HF dynamics. However, note that 
cotunneling current scales as I c t oc a 2 , whereas sequen- 
tial tunneling current I oc a; accordingly, cotunneling 
current can always be suppressed by making the tunnel 
barriers less transparentPH Moreover, inelastic cotun- 
neling processes exciting the QD spin can be ruled out 
for eV, ksT < u> due to energy conservation^. The 
effectiveness of a single quantum dot to act as an elec- 
trically tunable spin filter has also been demonstrated 
experimentally^ The spin-filter efficiency was measured 
to be nearly 100%, with I ct being smaller than the noise 
floor ~ 10 fA. Its actual value has been calculated as 
~ 10~ 4 fA, from which we roughly estimate I ct ~ 10~ 2 fA 
in our setting. This is smaller than the initial HF me- 
diated current I (t — 0) and considerably smaller than 
I mlXl even for an initially not fully polarized nuclear spin 
ensemble. Still, if one is to explore the regime where 
cotunneling cannot be neglected, phenomenological dis- 
sipative terms - effectively describing the correspond- 
ing spin-flip and pure dephasing mechanisms for inelastic 
and elastic processes respectively - should be added to 
Eqn.pH). 



VII. CONCLUSION AND OUTLOOK 

In summary, we have developed a master equation 
based theoretical framework for nuclear spin assisted 
transport through a QD. Due to the collective nature of 
the HF interaction, it incorporates intriguing many-body 
effects as well as feedback mechanisms between the elec- 
tron spin and nuclear spin dynamics. As a prominent 
application, we have shown that the current through a 
single electrically defined QD in the spin-blockade regime 
naturally exhibits superradiant behavior. This effect 
stems from the collective hyperfine interaction between 
the QD electron and the nuclear spin ensemble in the QD. 
Its most striking feature is a lifting of the spin-blockade 
and a sudden peak in the leakage current. The experi- 



mental observation of this effect would provide clear evi- 
dence of coherent HF dynamics of nuclear spin ensembles 
in QDs. 

Finally, we highlight possible directions of research go- 
ing beyond our present work: Apart from electronic su- 
perradiance, the setup proposed here is inherently well 
suited for other experimental applications like dynamic 
polarization of nuclear spins (DNP): In analogy to op- 
tical pumping, Eqn. ^ describes electronic pumping of 
the nuclear spins. Its steady states are eigenstates of A z , 
which lie in the kernel of the collective jump-operator A~ . 
In particular, for a completely inhomogeneous system the 
only steady state is the fully polarized one, the ideal ini- 
tial state required for the observation of SR effects. When 
starting from a completely unpolarized nuclear state, the 
uni-directionality of Eqn. ^ - electrons with one spin 
orientation exchange excitations with the nuclear spins, 
while electrons of opposite spin primarily do not - im- 
plies that the rather warm electronic reservoir can still 
extract entropy out of the nuclear system. More gener- 
ally, the transport setting studied here possibly opens up 
the route towards the (feedback-based) electronic prepa- 
ration of particular nuclear states in single QDs. This is 
in line with similar ideas previously de veloped in double 
QD settings, see e.g. Refs! 12 ' 15 ' 18 ' 2 " ' 3 ^. 

In this work we have specialized on a single QD. How- 
ever, our theory could be extended to a double QD 
(DQD) setting which is likely to offer even more possibil- 
ities. DQDs are routinely operated in the Pauli-blockade 
regime where despite the presence of an applied source- 
drain voltage the current through the device is blocked 
whenever the electron tunneling into the DQD has the 
same spin orientation as the one already present. The 
DQD parameters and the external magnetic field can be 
tuned such that the role of the states \a) , a =J,, f, in our 
model is played by a pair of singlet and triplet states, 
while all other states are off-resonant. Then, along the 
lines of our study, non-linearities appear due to depen- 
dencies between the electronic and nuclear subsystems 
and collective effects enter via the HF-mediated lifting of 
the spin-blockade. 

While we have focused on the Markovian regime and 
the precise conditions for the validity thereof, Eqn. ( 14 ) 



offers a starting point for studies of non-Markovian effects 
in the proposed transport setting. All terms appearing 
in the memory kernel of Eqn. ( 14 ) are quadratic in the 



fermionic creation and annihilation operators allowing for 
an efficient numerical simulation, without having to ex- 
plicitly invoke the flatness of the spectral density of the 
leads. This should then shed light on possibly abrupt 
changes in the QD transport properties due to feedback 
mechanism between the nuclear spin ensemble and the 
electron spin. 

Lastly, our work also opens the door towards studies 
of dissipative phase transitions in the transport setting: 
when combined with driving, the SR dynamics can lead 
to a variety of strong-correlatio n effe cts, non-equilibrium 
and dissipative phase transition d 1 * 57 * ^, which could now 
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be studied in a mesoscopic solid state system, comple- 
menting other approaches to dissipative phase transitions 
in QDsEsHSl. 
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Appendix A: Microscopic Derivation of the Master Equation 

In this Appendix we provide some details regarding the derivation of the master equations as stated in Eqn.Q 
and Eqn.(28|. It comprises the effect of the HF dynamics in the memory kernel of Eqn.(12| and the subsequent 



approximation of independent rates of variation. 

In the following, we will show that it is self-consistent to neglect the effect of the HF dynamics L\ (t) in the memory- 
kernel of Eqn.(12| provided that the bath correlation time r c is short compared to the Rabi flips produced by the 



HF dynamics. This needs to be addressed as cooperative effects potentially drive the system from a weakly coupled 
into a strongly coupled regime. First, we reiterate the Schwinger-Dyson identity in Eqn.(13l as an infinite sum over 
time-ordered nested commutators 
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More explicitly, up to second order Eqn.(|Al|) is equivalent to 
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Note that the time-dependence of Hi (r) is simply given by 



Hi(t) 



-H_+H AOR , H ± = 9 -S ± AT 1 



(A3) 



(A4) 



where the effective Zeeman splitting uj = uj + g (A z ) t is time-dependent. Accordingly, we define C\ (t) = C + (r) + 
£_ (r) + £aoh (t) = e luJT C + + e~ lUT C- + £aoh, where C x - — [H x , •] for x = ±, AOH. In the next steps, we will 
explicitly evaluate the first two contributions to the memory kernel that go beyond n = and then generalize our 
findings to any order n of the Schwinger-Dyson series. 



First order correction 



The first order contribution n = 1 in Eqn.(12| is given by 

ft 



S« =ij dr j driTr B (c T e- lC ° T [ffi (n) , a]) 
Performing the integration in n leads to 
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where, for notational 

e -iHor [ Hr ps (t-T)p%\ e m ° 



re introduced 
H T (r),ps(t)p° B . 



r)p% 



and X T 



convenience. we unnxinced the operators X — CtPS (t 

In accordance with previous approximations, we have re- 

place d c m,T p s (t - r) e lH ° T by p s (t) since any additional term besides H would be of higher order in perturbation 
theorySIM. j n particular, this disregards dissipative effects: In our case, this approximation is valid self-consistently 
provided that the tunneling rates are small compared to effective Zeeman splitting u). The integrand decays on the 
leads-correlation timescale r c which is typically much faster than the timescale set by the effective Zeeman splitting, 
lot c <C 1. This separation of timescales allows for an expansion in the small parameter lot, e.g. - {e lulT — 1J « igr. 
We see that the first order correction can be neglected if the the bath correlation time r c is sufficiently short compared 
to the timescale of the HF dynamics, that is gr c <C 1. The latter is bounded by the total hyperfine coupling constant 
Ahf (since Hg-A^H < Ahf) so that the requirement for disregarding the first order term reads Ahft c <C 1. 



Second order correction 



The contribution of the second term n — 2 in the Schwinger-Dyson expansion can be decomposed into 

=(2) _ n(2) 
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The first term contains contributions from Haoh only 
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Similarly, ^ s which comprises contributions from Hg only is found to be 
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Here, we have used the following simplification: The time-ordered products which include flip-flop terms only can be 
simplified to two possible sequences in which £ + is followed by £_ and vice versa. This holds since 



C ± C ± X = [H±, [H±,X]] = H ± H ± X + XH ± H ± - 2H ± XH ± = 0. 



(All) 



Here, the first two terms drop out immediately since the electronic jump-operators S fulfill the relation S ± S ± = 0. 



In the problem at hand, also the last term gives zero because of particle number superselection rules: In Eqn.( 12 1 the 
time-ordered product of superoperators acts on X = [Ht, ps {t ~ t) p%\ . Thus, for the term H±XH± to be nonzero, 
coherences in Fock space would be required which are consistently neglected; compare RefP^. This is equivalent 
to ignoring coherences between the system and the leads. Note that the same argument holds for any combination 
H^XHv with /i, v = ±. 

Similar results can be obtained for Sf ' which comprises H± as well as i?AOH in all possible orderings. Again, using 
that the integrand decays on a timescale r c and expanding in the small parameter lot shows that the second order 
contribution scales as ~ (gT c ) 2 ■ Our findings for the first and second order correction suggest that the n-th order 
correction scales as ~ (gT c ) n . This will be proven in the following by induction. 



n-th order correction 



The scaling of the n-th term in the Dyson series is governed by the quantities of the form 
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where the index suggests the order in which H± (giving an exponential factor) and 7?aoh (resulting in a factor 
of 1) appear. Led by our findings for n — 1,2, we claim that the expansion of ^1 (r) for small lot scales as 
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( r ) ~ (s T )™- Then, the (n + l)-th terms scale as 
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Since we have already verified this result for n = 1,2, the general result follows by induction. This completes the 
proof. 



Appendix B: Adiabatic Elimination of the QD Electron 



For a sufficiently small relative coupling strength e the nuclear dynamics are slow compared to the electronic QD 
dynamics. This allows for an adiabatic elimination of the electronic degrees of freedom yielding an effective master 
equation for the nuclear spins of the QD. 

Our analysis starts out from Eqn.(33l which we write as 



p = W p + Wip, 



where 



W p = -i[oj a S z ,p] +7 
Wip = -i[H HF ,p]. 



s- P s + --{s + s-,p} 



S Z P S Z - -p 
4 



(Bl) 

(B2) 
(B3) 

Note that the superoperator Wo only acts on the electronic degrees of freedom. It describes an electron in an external 
magnetic field that experiences a decay as well as a pure dephasing mechanism. In zeroth order of the coupling 
parameter e the electronic and nuclear dynamics of the QD are decoupled and SR effects cannot be expected. These 
are contained in the interaction term Wi. 

Formally, the adiabatic elimination of the electronic degrees of freedom can be achieved as follows^: To zeroth 
order in e the eigenvectors of Wo with zero eigenvector Aq = are 



WoA* ® Pss = 0, 



(B4) 



where pss = 14-) (44 IS the stationary solution for the electronic dynamics and p describes some arbitrary state of 
the nuclear system. The zero-order Liouville eigenstates corresponding to Ao = are coupled to the subspaces of 
"excited" nonzero (complex) eigenvalues A& ^ of Wo by the action of Wi . Physically, this corresponds to a coupling 
between electronic and nuclear degrees of freedom. In the limit where the HF dynamics are slow compared to the 
electronic frequencies, i.e. the Zeeman splitting loq, the decay rate 7 and the dephasing rate T, the coupling between 
these blocks of eigenvalues and Liouville subspaces of Wo is weak justifying a perturbative treatment. This motivates 
the definition of a projection operator P onto the subspace with zero eigenvalue Ao = of Wo according to 



Pp = Tr e , [p] <g) pss = H O II) (II 



(B5) 



where p = Tr e \ [p] is a density operator for the nuclear spins, Tr e | . . . denotes the trace over the electronic subspace 
and by definition WqPss — 0. The complement of P is Q = 1 — P. By projecting the master equation on the P 
subspace and tracing over the electronic degrees of freedom we obtain an effective master equation for the nuclear 
spins in second order perturbation theory 

p = Tr e , [PWxPp - PWiQW^QWiPp] ■ (B6) 

Using Tr e | [S z pss] = — 1/2, the first term is readily evaluated and yields the Knight shift seen by the nuclear spins 



Tr el [PWiPp] =+i y -[A z ,p). 



The derivation of the second term is more involved. It can be rewritten as 

- Tr e , [PWxQW^QWxPp] = -Tr el [PWi (1 - P) W^ 1 (1 - P) WiPp] 

poo poo 

= dTTr e | [PWie WaT WiP P ] - / dr Tr e , [PWiPWiPp] . 
Jo Jo 



(B7) 

(B8) 
(B9) 
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Here, we used the Laplace transform — W 1 = f n ° ° dr e w ° T and the property e WoT P = Pe WoT = P. 
Let us first focus on the first term in Eqn.( |B9| . It contains terms of the form 

Tr el [P [A+S-,e WoT [A~S + , p ® p ss ]]] = Tr el [S- e w ° T (S + Pss )] A+A~p 

-Tr el [S-e w ° T (S+pss)] A~ pA+ 
+Tr el [S-e w ° T ( PSS S + )] pA~ A+ 
-Tr el [S-e w ° T ( PSS S+)] A+ pA~ 

This can be simplified using the following relations: Since pss = \i) (41 j we have S~pss = and pssS + = 0. 
Moreover, |f) (4| and |4) (f| are eigenvectors of Wo with eigenvalue — (iuj + a/2) and + (iuj + a/2), where a = 7 + L, 
yielding 



(BIO) 
(Bll) 
(B12) 
(B13) 



This leads to 



e WoT (S+Pss) = e-^ +a ' 2 ^ |t> (II 
e WoT (pssS-) = e +(-o+«/2)r j;) (t | 



Tr e , [P [A+S^e 14 ^ [A-S^,|t®/» ss ]]] = e"^ 0+Q / 2 ) r (A+A'p - A~ pA+) 



Similarly, one finds 



Tr e , [P [A~S + ,e WoT [A+ S~ , p ® Ass]]] = e + ^ 0+Q / 2)r (/x^+A" - . 



(B14) 
(B15) 

(B16) 
(B17) 



Analogously, one can show that terms containing two flip or two flop terms give zero. The same holds for mixed 
terms that comprise one flip-flop and one Overhauser term with ~ A Z S Z . The term consisting of two Overhauser 
contributions gives 



Tr e , [P [A z S z ,e w ° T [A Z S Z ,p ® pss]]] = -\ [2A Z pA z - [A z A z ,p]} 



(B18) 



However, this term exactly cancels with the second term from Eqn.(B9|. Thus we are left with the contributions 
coming from Eqn.(B16) and Eqn.(B17|. Restoring the prefactors of —ig/2, we obtain 

„2 



Tr el [PWiQ (-Wq" 1 ) QWiPp] = 



-(iuj +a/2)T 



(A-pA+ - A+A-p) 



(A-pA+ - pA+A- 



Performing the integration and separating real from imaginary terms yields 



(B19) 



Tr el [PW^ (-n^ 1 ) QWiPp] 



A-pA+ - l -{A+A-,p} 



+ ia [A+A , p] , 



(B20) 



where c r = g 2 / (4w 2 + a 2 ) a and a = g 2 / (4w 2 + a 2 ) uj a . Combining Eqn.(B7) with Eqn.(B20) directly gives the 
effective master equation for the nuclear spins given in Eqn.([3]) in the main text. 
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